A proposal to determine the spectrum of pairing-glue in high-temperature 

superconductors 
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We propose a method for an analysis of the angle-resolved photoemission data in two-dimensional 
anisotropic superconductors which directly yields the spectral function of the bosons mediating 
Cooper pairing. The method includes a self-consistency check for the validity of the approximations 
made in the analysis. We explicitly describe the experimental data needed for implementing the 
proposed procedure. 



Introduction. Understanding the mechanism of su- 
perconductivity in the high temperature superconducting 
cuprates requires knowledge of the spectral function of 
the bosons mediating the pairing as well as their coupling 
to the fermions. For conventional superconductors, such 
as Pb or Al, the phonon mechanism of superconductivity 
was firmly established through development of a proce- 
dure to invert the tunneling spectra (conductance as a 
function of voltage) and obtain the phonon spectral den- 
sity 0, 0] • Here we propose a generalization of this pro- 
cedure for anisotropic superconductors, which in general 
may have pairing mediated by the electron-electron inter- 
actions, using the angle resolved photoemission (ARPES) 
data as an input. 

Recall that in conventional superconductors the de- 
pendence of the normalized tunneling conductance, 
G s (V)/G n (V), where s and n refer to the superconduct- 
ing and the normal state respectively, on the energy de- 
pendent gap function, A(u>), is known. In turn, A(w) 
is related to the spectrum of pairing bosons through the 
Eliashberg theory This makes the inversion of the 
tunneling spectra possible. 

Several important differences arise in the case of uncon- 
ventional superconductors. First, since for an anisotropic 
gap function the tunneling current strongly depends on 
direction, and since the tunneling spectra are integrals 
over all the directions (often with a varying weight due to 
the spatial structure of the matrix element), such spec- 
tra in general neither provide sufficient information to 
resolve the anisotropy of the boson-fermion coupling nor 
allow its frequency dependence to be determined quan- 
titatively. Second, in contrast to the case of phonons, 
in superconductors where the pairing is mediated by 
electron-electron interaction, the pairing spectra change 
and are temperature dependent below the transition tem- 
perature, T c , when the single particle spectra are gapped. 
Also, since the typical energy of interaction between elec- 
trons is high, any inversion procedure must determine 
whether Eliashberg theory is suitable for the analysis. 

We propose an inversion procedure using the ARPES 
data, which is applicable to anisotropic superconductors, 
uses the Eliashberg theory and contains a consistency 
check 0. Below, we first summarize the procedure and 



the experimental input required to implement it, and 
then describe the technical details. 

Outline of the procedure. As input the following ex- 
perimental information for an optimally doped or over- 
doped sample is required: (i) For Normalization. Just 
above T c , for several points at the Fermi surface, we re- 
quire the integrated area under the momentum distribu- 
tion curve (MDC) for momentum perpendicular to the 
Fermi surface, as defined in Eq. (|5j) below. This step gives 
the anisotropic equivalent of normal state conductance in 
the McMiUan-Rowell procedure and provides the neces- 
sary normalization for the second step, (ii) For Inver- 
sion. With the same polarization and incoming photon 
energy, and at the same set of points at the Fermi sur- 
face, we need, at several temperatures below T c , the dif- 
ference in the intensity of the MDCs between the normal 
state above T c and the superconducting state, for ener- 
gies up to several times the superconducting gap. These 
data may be cither a function of momentum or (better 
still) integrated during measurement over the momentum 
normal to the Fermi surface as in (i). This step allows 
the extraction of the spectral density of the boson and 
the form of the coupling function. The data are used 
to extract the superconducting gap function, according 
to Ea. i|10[l . The gap function is then inverted by using 
Eliashberg theory to obtain the boson spectrum and its 
coupling to fermions in different angular channels accord- 
ing to Eq. (|18|) . The self-consistency check on the Eliash- 
berg theory is obtained by comparing the energy range 
of the bosons with the input ARPES spectra and with 
the estimates of the bandwidth. 

ARPES spectra. The intensity of the ARPES signal, 
7(k, to), for a fixed momentum and binding energy of the 
emitted electron, is 



/(k,w)=A<(k)/(«/T)A(k > a;) 



(1) 



Here / is the Fermi function, and M. is a prefactor, which 
depends on the square of the dipole matrix element be- 
tween the initial and the final states, and on the kinemat- 
ical factors. The information about the fermion spectrum 
is contained in the spectral function 



A(k,u) = --Jm G R (k,oj), 

7T 



(2) 
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where G R is a retarded Green's function. It is the k- 
dependent prefactor, M., that obscures the true momen- 
tum dependence of the spectral function, and prevents 
direct inversion of the ARPES intensity. 

In disentangling the momentum dependence of the 
spectral density, A(k, u>), from the signal intensity mea- 
sured in ARPES, we rely on several properties of the pref- 
actor. First, while M. depends on the momentum, k, it 
does not contain a significant dependence on the temper- 
ature, T, or the electron binding energy, uj Second, 
since the momentum dependence of M. (k) is determined 
by the spatial structure of the electron wave functions, 
the energy and momentum range over which it varies are 
of the order of the Fermi energy and the Brillouin Zone 
respectively. Consequently, if, for fixed iv, the spectral 
density, A, is sharply peaked in the momentum space, 
the dependence of the prefactor on |k| may be ignored. 

Step (i): MDCs and normalization. We write the mo- 
mentum k = (kf + fcj_, s), where fc/(s) is the Fermi mo- 
mentum and s is the tangential component of k at a point 
on the Fermi surface. The normal state Green's function 
at Matsubara frequencies, u> n — 2irT(n + 1/2), is 



Gjv(w n ,k) = [iu n Z(k,uj n ) - £(k)]~ 



(3) 



where Z contains the effects of the interactions, £(k) ss 
Vf(s)k±, and Vf is the bare, unrenormalized, Fermi ve- 
locity. Therefore the spectral function is 



A N (k,oj) = - 



w2"(k,w) 



ir (wZ'(k, u) - v f k ± ) 2 + (tuZ"(k, w)) 2 ' 

(4) 

where Z = Z' + iZ" is the analytic continuation of 
Z(k 7 uj n ) to real frequencies. 

Our subsequent analysis assumes that Migdal's theo- 
rem is valid for high-T c superconductors, and checks for 
the validity of this assumption at the end. The essential 
statement of the theorem is that, when the energy range 
of the interaction is small compared to the electron band- 
width but the momentum ranges of the electrons and 
the boson are comparable, (a) the electronic self-energy 
is momentum-independent, (b) the vertex corrections to 
the interactions are negligible 6]. Under such conditions, 
we can approximate Z(k,u> n ) w Z(s,uj n ; kf) in Eq.Q to 
a high degree of accuracy. We believe that this assump- 
tion will be borne out by the self-consistency check since 
both the neutron scattering data and the analysis of 
the ARPES data in the normal state 0, Q indicate the 
range of the interactions to be 0.2-0.4 eV, only a fraction 
of the bandwidth. 

In contrast to the usual Eliashberg theory, we must 
consider the change in the spectrum of the boson in the 
superconducting state. We show below that the inversion 
of the ARPES data at different temperatures below T c 
explicitly gives the evolution of the pairing interaction as 
the superconducting order develops. This point is very 
important because the value of T c is determined by the 



spectrum just above T c , while all the properties below 
T c (for example the temperature dependence of the su- 
perconducting gap) are determined by the self-consistent 
change of the pairing spectrum as a function of temper- 
ature. 

The momentum distribution curves (MDC), are ob- 
tained by plotting the intensity of the ARPES spectra, 
J(k, u>) at a fixed electron energy to as a function of k 
|8|. The MDCs needed here are those taken along the 
lines normal to the Fermi surface for several direction. 
These curves have Lorentzian shape as a function of k± 
H, which confirms (by inspection of Eq.@) that both 
the renormalization function, Z, and the prefactor M 
depend only weakly on fcj_ in the neighborhood of the 
peak of the MDC, M(s,k f + k±) « M(s, kf) = M(s). 

m 

Divide the ARPES intensity by the Fermi function to 
remove the temperature dependence, and consider the 
area under thus rescaled MDC, 



dk 



7(«/T) 



M(s) / dk ± A(k,u). (5) 



It follows from the Lorentzian shape of the spectral func- 
tion in the normal state that 



Jn(s,lo) 



Vf(s) 



(0) 



is independent of the binding energy, and only depends 
on the position on the Fermi surface via the prefactor 
and the bare Fermi velocity. Therefore it serves to nor- 
malize the MDCs. It also follows that the energy and 
temperature dependence of the area under the MDC is 
simply given by the Fermi function. Note that to recover 
most of the area under the MDCs in the experiment it is 
sufficient to carry out the integration over the range of 
k± of only a fraction of the Brillouin Zone. 

Step (ii): Extraction of the gap function. Next con- 
sider the difference between the electron Green's func- 
tions in the normal and the superconducting state, and 
employ the spectral representation to write this difference 
as 

Gjv(k,w n ) - G s (k,w n ) (7) 
" +0 ° dx 



(A N {k, x)-A s (k,x)y 



Here G s is the '11' component of the matrix (Nambu- 
Gorkov) Green's function of a superconductor 



G(k,cj„) = - 



iu n Z(k,Lj n ) +£(k)r 3 + $(k,^„)(ir 2 ) 

Z 2 (k,LO n )L0 2 n + ^(KLO n )+e 



(8) 

and <& is the anomalous part of the self-energy. Wc 
now integrate both parts of Eq. Q over the momentum 
normal to the Fermi surface at a given point s. Using 
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Migdal's theorem, we replace both $ and Z in the inte- 
grand by their values at the Fermi surface, and introduce 
the gap function A(s,u> n ) = ^(s,u) n )/Z(s,u> n ) to find 



ITT 



sgn(w„) 



U)„ 



f +oa dx 



(9) 



J dk±(A N (k,x) -A s (k,x] 



Assuming particle-hole symmetry, A(—kj_,s,—x) = 
A(k±, s, x), and the existence of the center of inversion 
for the Fermi surface, Ai(k) = M.{— k), we rewrite the 
integral over the energy, x, via an integral over only the 
occupied states, x < (cf.0]). Using the normal state 
result, Eq.©, to eliminate M(s)/vf(s), we find (oj n > 0) 



A(s,u n )=u) n [l-w n K(u n ,s)] 2 -l 



1/2 



(10) 



2 i 

K.(u n ,s) = - [Jn(s)] 



(11) 



dx 



UJ 2 



\Jn{s,x) - J a (a,x)] 



At each uj n = irT(2n + 1) and for given direction s at 
the Fermi surface, the function K(u) n , s) depends solely 
on experimentally measured intensities. Therefore us- 
ing Eas. i|lUfl - Hllfl the gap function A(s ) u> n ) can be de- 
termined directly from the ARPES spectra. 

Inversion of the gap equation. We assume a pairing 
function of the form g(k, k')| 2 _B(k — k', u> n — w n '), where 
g(k, k')| is the matrix element of the interaction, and 
.B(k, u) n ) is the propagator of the boson. In the Eliash- 
berg theory the equations for the normal and the anoma- 
lous part of the self-energy take the form 



i -z(k, w „)W„ = -tJ2 dk'| g (k,k') 



ioy Z(k', ay ) B (k - k' , uj n - u) n > ) 
u>l, Z 2 (k', Un , ) + $ 2 (k', ^ ) + C 2 (k' ) 

$ 2 (k', w „0+e 2 (k')' 



(12) 
(13) 



where 77 = +1 (77 = — 1) for coupling in the channel even (odd) under time-reversal |l2|. Integrating over k' ± in the 
right hand side, we arrive at an equation for the gap function 



uj n A(s,u} n ) = -7rT^ J ds'g(s,s')B(s - s',u n - w n i)[G(s,w n ,)A(s,w n ) -riw n F(s,u n >)], 



(14) 



r 



where g(s,s') = |g(s, s')\ 2 /\vf(s)\ is the effective coupling 
at the Fermi surface, and G and F are the ^-integrated di- 
agonal and off-diagonal components of the Green's func- 
tion respectively, 



G(s,uj n ) = 
F(s,uj n ) = 



+ A 2 (s,uj n ) 
A(s,uj n ) 
+ A 2 (s,w n )' 



(15) 
(16) 



We now make a simplifying assumption (to be verified) 
that the pairing function is separable and diagonal in the 
space of the Fermi surface harmonics, ipi(s). These har- 
monics are the orthonormal eigcnfunctions of the sym- 
metry operators of the Fermi surface, and correspond to 
different angular channels for a cylindrical Fermi surface 
Therefore 

g(s, s')B(s - s', uj n - uv) (17) 



Expand in ipi(s) all functions of s in Ea. (|14|l . A(s,u) n ) = 
J2i Ai(w„)^i(s) and the same for functions G and F, to 
obtain a linear system for the coupling function 



{3ik(Un,Un - Wn')-Bfe(uv) = OCi(u) n ). 



(18) 



All the coefficients in this equation are known from the 
form of the gap function and the Fermi surface, 



f3 lk (uj n ,uj n ,) = irT T]6 lk uj n F k (uj n ,) - b ik {uj n )G k {uj n ,) 
hk{u n ) = y^aijfeAj(o;„), aj(w„) = w n Ai(u) n ), 



dsipi(s)ijjj{s)ijjk(s). 



Solution of Ea. (|18|l gives the pairing function at Mat- 
subara frequencies and in particular angular channels. 
It can then be analytically continued to real frequencies 
using either the Pade approximation or by solving 
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the integral equation obtained by analytic continuation 
of Ea. l(TH|l in analogy to Reflla. The full pairing function, 
J2 i Bi(u))if>i(s)il}i(s') can then also be used to obtain the 
self-energies, Z and using Eas.lfT2 |l -l(T3 |l . 

Self- consistency. Two major assumptions have been 
made in our analysis, and have to be checked for self- 
consistency. First, we assumed that Migdal's theorem 
holds. If the energy range of the pairing function is sig- 
nificantly smaller than the bandwidth in all channels i, 
the procedure is self-consistent. Second, we assumed sep- 
arable form of the pairing function, Ea. ll7|) . In a realistic 
calculation we truncate the number of angular channels, 
i, in Ea. (|18(l : if inclusion of an additional channel does 
not change appreciably the pairing functions in the other 
channels, the approximation is justified. Note that a min- 
imum of two channels is required: the isotropic chan- 
nel which alone determines the self-energy in the normal 
state if the Fermi-surface is circular, and the d x i_ y i chan- 
nel responsible for pairing; we expect additional channels 
to enter with smaller couplings. The change of spectral 
strength with temperatures in these channels, obtained 
via the outlined procedure, yields crucial information on 
the change of self-energies below T c . 

Experimental constraints. Most available ARPES 
spectra are taken on B^S^CaC^Os+s samples. Since 
our analysis includes only superconducting order and 
does not account for the pseudogap, the measurements 
should be done on an optimally doped or overdoped sam- 
ple. In the normal state the width of the peak in MDCs 
is a fraction of the Brillouin Zone || , and we expect that 
the data in the superconducting state have to be collected 
over the same window in momentum space. Since the 
area Jn in Eq.|0 does not depend on the energy to, and 
since MDC is a Lorentzian with a flat (^-independent) 
background at to = 0, we suggest that Jn is most easily 
determined at the Fermi energy. Away from the Fermi 
energy the background becomes momentum-dependent; 
however, the difference between the signal in the normal 
and superconducting state, required in Eq.JHJ is insen- 
sitive to the background. In that equation we expect 
the energy range over which the k± -integrated difference 
in intensity is appreciable to be of the order of a few 
times the superconducting gap, so that the data should 
be taken up to to ~ 0.1 — 0.2eV. 

Such data may be difficult to take near the (ir, 0) points 
because of the bilayer splitting, not accounted for in our 
one-band approach (generalization to a multi-band case 
is, however, straightforward). Also, near these points 
the MDCs corresponding to different Fermi arcs begin to 
overlap at energies away from to = 8J. Consequently, 
it is more feasible to measure the spectra close to the 
nodal directions; we expect that spectra taken over half 
of the angular range between the node and the antinode 
would contain enough information to carry out the anal- 
ysis suggested here. Since the analysis of the spectra is 
complicated by the superstructure in Bi-0 layer, sam- 



ples containing Pb, which suppresses Bi-0 modulations, 
may be more suitable. Finally, the resolution function of 
the measurement apparatus has to be disentangled from 
the signal intensity; the separation can be done by an 
integral transform if the energy dependence of the res- 
olution function is known. A numerical procedure for 
implementing the method has been developed and will 
be made available. 

We have proposed a method to extract the spectrum of 
the pairing function from the ARPES measurements in 
high temperature superconductors which makes a mini- 
mum of theoretical assumptions, does not use any prej- 
udicial form of the pairing fluctuation spectrum as an 
input, and contains a self-consistency test. In contrast 
to the analysis of optical conductivity , our approach 
uses single particle spectra which contains more infor- 
mation and explicitly gives the evolution of the pairing 
spectrum with temperature in the superconducting state. 
Implementation of this approach through appropriate ex- 
periments promises to yield the answer to one of the cen- 
tral questions in the high T c problem. 
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